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Abstract 

We study the equilibrium properties of a lattice-gas model of an A + B — > catalytic 
reaction on a one-dimensional chain in contact with a reservoir for the particles. The 
particles of species A and B are in thermal contact with their vapor phases acting as 
reservoirs, i.e., they may adsorb onto empty lattice sites and may desorb from the lattice. 
If adsorbed A and B particles appear at neighboring lattice sites they instantaneously 
react and both desorb. For this model of a catalytic reaction in the adsorption-controlled 
limit, we derive analytically the expression of the pressure and present exact results for 
the mean densities of particles and for the compressibilities of the adsorbate as function 
of the chemical potentials of the two species. 



PACS: 05.50.+q; 64.60.Cn; 68.43.De; 82.65.+r 



1 Introduction 



During the two decades following the work of Ziff, Gulari, and Barshad (ZGB) [1] there has 
been a remarkable development in the theoretical analysis of catalytically activated reactions. 
The ZGB model, sometimes referred to as the "monomer-dimer" model, has been introduced 
to describe the important process of oxidation of carbon monoxide on a catalytic surface [2]. 
Within this model, a monomer (CO) adsorbs onto a single vacant site of the surface (since 
dissociation of CO is not considered, CO is treated as a monomer), while a dimer (O2) adsorbs 
onto a pair of adjacent vacant sites and then immediately dissociates. Both CO and O2 are in 
thermal contact with their gaseous reservoirs. Nearest neighbors of adsorbates composed of 
a dissociated O atom and a CO molecule react and form a CO2 molecule, which then desorbs 
from the metal surface. The ZGB model predicts remarkable new features compared to the 
classical mean-field kinetic schemes of rate equations [2]: for a two-dimensional catalytic 
substrate as the CO gas pressure is lowered the system undergoes a first-order transition 
from a CO saturated inactive phase (zero rate of CO2 production) into a reactive steady state 
(non-zero rate of CO2 production) followed by a continuous transition into an 02-saturated 
inactive phase. This continuous transition was shown to belong to the same universality class 
as the directed percolation and the Reggeon field theory [3]. A simpler "monomer-monomer" 
model, in which particles of both species can adsorb on single but different sites, has also been 
proposed [4-7], and it has been shown that in (2+1) spatial dimensions this system exhibits a 
first-order transition from a phase saturated with one species to one saturated with the other. 
Allowing desorption of one species leads to a continuous transition which also belongs to the 
directed percolation universality class [7, 8] . 

These observations have significantly increased the interest in the properties of such models 
of catalytically activated reactions. Different aspects of the dynamics of the adsorbed phase 
for these two models, as well as for their extensions to molecules with more complicated 
structures (e.g., dimer-dimer [9], dimer-trimer [10] models), have been examined thoroughly 
[11-23], and the studies have confirmed an essentially collective, many-particle behavior. On 
the other hand, the equilibrium properties of catalytically activated reactions in systems 
in which the reactive species undergo continuous exchanges with their vapor phases acting 
as reservoirs, i.e., adsorb on and desorb from the catalytic substrate, have been much less 
studied, and, consequently, the understanding of the equilibrium state remains rather limited. 
Only recently such equilibrium properties, revealing a rather non-trivial behavior, have been 
obtained for the more simple case of single-species reactions A + i^0ona one-dimensional 
chain with a random distribution of catalytic segments or catalytic sites [24,25]. Kinetics of 
the diffusion-limited A + A — > reactions, controlled by the constraint that the particles may 
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undergo reactions only when they meet each other in the vicinity of special catalytic sites has 
been discussed in Ref. [26]. 

Here we study the properties of the equilibrium state for a very simple model of a 
monomer-monomer A + B — > catalytic reaction on a one-dimensional chain. Within this 
model the particles of species A and B undergo continuous exchanges with the mixed vapor 
phase, i.e., they adsorb onto empty lattice sites on the chain and may thermally desorb from 
them, while the vapor phase as a reservoir is steadily maintained at constant chemical poten- 
tials. In addition, if any adsorbed A and B particles appear at neighboring lattice sites, they 
instantaneously react and both leave the chain (i.e., desorb), and the reaction product AB is 
completely removed from the system. For this model we derive the expression for the pressure 
of the adsorbed particles, which provides the complete thermodynamical description of the 
system, and we present exact asymptotic (in respect to the chemical potential) results for the 
mean densities of the particles and for the compressibilities of the system. We note that one 
could also consider a model in which the product AB is not removed from the system but 
is forming its own vapor phase from which it may readsorb on the substrate. However, this 
more complicated case will not be studied here. 

In Sec. 2 we describe the model and introduce basic notations. In Sect. 3 we derive the 
recursion relations obeyed by the partition function of the system and obtain the solution of 
these recursion relations. The connection of this model with classical spin S = 1 models is 
discussed in Sect. 4. In Sect. 5 we derive the expressions of the mean densities of the particles 
and of the compressibilities of the system, analyze their asymptotic behavior as function of 
the chemical potentials, and present results of Monte Carlo simulations for the occupation 
of the lattice in the limit of infinite chemical potentials. We conclude the paper with a brief 
summary in Sect. 6. 

2 Model 

Consider a one-dimensional, regular lattice of N adsorption sites (Fig. 1), which is brought 
in contact with the mixed vapor phase of two types A and B of hard-core particles without 
attractive interactions. The difference in chemical potential between the vapor phase and 
the adsorbed phase ha and fi B , including the binding energy of an occupied site, is main- 
tained constant (thus HA,B > corresponds to a preference for adsorption), and the constant 
activities are hence defined by 

z A = exp((3fi A ), z B = exp((3(j, B ), (1) 

where = k B T is the thermal energy, T is the temperature, and k B is the Boltzmann 
constant. 
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Figure 1: One-dimensional regular lattice of N adsorption sites in contact with a mixed vapor phase. 
Black and gray circles denote A and B particles, respectively, (a) denotes a configuration in which an 
instantaneous reaction takes place upon which both particles desorb and the product AB leaves the 
system. The sites i = and i = N + 1 are always empty. 



The A and B particles can adsorb from the vapor phases onto vacant sites, and can desorb 
back to the reservoir. The occupation of the i-th. site is described by two Boolean variables 
rii and mj, such that 

j 1, if the i-th. site is occupied by an A particle, 

Ui = \ 

y 0, otherwise, 

and 

I 1, if the i-th site is occupied by a B particle, 
rrii = < 

I 0, otherwise, 

and thus the state of the site i is specified by (ni,mi) £ {(0, 0), (0, 1), (1, 0), (1, 1)}. Note 
that the state (1, 1), which in general is allowed, has a zero probability in the case of hard- 
core interaction. (The equivalent description in terms of <Tj = 0, +1,— 1 for empty sites and 
sites occupied by A and B, respectively, will be discussed in Sect. 4.) For computational 
convenience we add two extra "boundary" sites i = and i = N + 1 at both ends of the chain 
and stipulate that these sites are always unoccupied, i.e., no = njy+i = ra^ = uin+i = 0. 
Furthermore, the system is locally "frustrated", i.e., the equilibrium between the chain and 
the reservoir is locally perturbed by the "kinetic" constraint that as soon as an A and a B 
particle occupy neighboring sites they instantaneously react and leave the chain (and the AB 
product leaves the system); therefore, our model corresponds to a purely adsorption controlled 
dynamics [15]. 

In equilibrium the partition function of this two-species adsorbate on a chain of N sites, 
in which the configurations {(ni,nii)} with neighboring A and B or double occupancy of a 
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site are excluded, can be written as: 

N 

z »= e n#*B**' (2) 

{(n„mi)} i=1 

where <3?j is an indicator function with the following properties: 
1, if site i is empty, 

1, if site i is occupied by an A particle, i.e., (rij,mj) = (1,0), while rrii±i = 0, 
1, if site i is occupied by a B particle, i.e., (n^, mj) = (0, 1), while n^±i = 0, 
0, otherwise. 

A very simple realization of the indicator function <£j is given by 

$i = (l - riim^j (l - nirrii-ij (l - riim i+1 ^ (l - m^rai-i) - m^nj+i) . (3) 

The factor (1 — niirii) ensures that each site i can be occupied by either A or B particle, or 
be empty, while the product ^1 — niirii-i^ (l — nj?7ij + i^ ensures that if the site i is occupied 
by an A the neighboring sites are not occupied by a B (and vice versa for the last two terms). 
Note that <E>j in Eq. (3) is symmetric with respect to the variables rij and mj, and therefore 
the partition function in Eq. (2) is a symmetric function of za and zb- 

The thermodynamics is given by the pressure (in units of the square of the lattice spacing) , 

P = P{T,vL A ,vL B ) = ±.Vm (4) 

p N— >oo iv 

from which all other quantities of interest, such as the mean densities of the species A and B 
and the compressibilities can be obtained straightforwardly by differentiating it suitably with 
respect to the chemical potentials [ia and \i B - 

Before closing this section we mention two trivial cases: 

(a) If the A and B particles do not react and do not interact, such that a lattice site can be 
occupied by an i, a S, or simultaneously by an A and a B particle, i.e., <E>j = 1, the pressure 
and the mean density of the A (B) particles is given by the classical Langmuir adsorption 
model 

P =- In (1 + z a ,b) and n A ,B = : ZA ' B ■ (5) 
P 1 + ZA,B 

(b) If the A and B particles do not react, but obey the hard-core exclusion interaction 
such that a lattice site cannot accommodate more than one particle (either A or B), i.e., 
$i = (1 — niirii), one has 

P = \ In (1 + z A + z B ) and n A ,B = : , ZA,B , • (6) 

(3 1 + ZA + z B 

For the following these results will serve as reference expressions, allowing us to identify that 
contribution to the thermodynamics of the system which stems from the "reaction" part 
alone. 
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3 Calculation of the partition function based on recursion re- 
lations 



We define two auxiliary partition functions: 



7 (A) 



N 



{(ra^mf)} i=l 



and Z™ = U z A z b^ 



N 



n N =l 



{(rii,mi)} i=l 



(7) 



mjv=l 



which are the partition functions of the chain under the constraint that the site i = N 
is occupied by an A particle (i.e., njy = 1) or a B particle (i.e., tun = 1), respectively. 
Summing over the occupation variables of the site i = N yields the following relation for the 
unconstrained partition function (Eq. (2)): 



Z N = Z N _ X + + Z { *\ for N > 2. 

Summing over the occupation variables of the site i = N — 1 for njv = 1 or mj/ 
respectively, one finds that the constrained partition functions (Eq. (7)) satisfy 



7(C) 

J N 



z c Z N _ 2 + z c Z^\, for N > 3, C = A, B. 



(8) 
■ 1, 

(9) 



Equations (8) and (9), augmented with the obvious initial conditions 



Zi = 1 + z A + z B , 



MA) _ ry(A) _ . 2 

Zj 1 — ZA, ^2 — Z A + %Ai 



7(B) _ 



Z B, 



7(B) _ 



ZB + Z B , 



(10) 



represent a closed system of coupled recursion relations determining Z^ for arbitrary N , za, 
and z B . 

The solution of Eqs. (8)- (10) can be found using the standard generating function tech- 
nique. Let 



z 



(C) 



Y zPt N , C = A,B, 



(11) 

N=l N=l 
denote such generating functions. Multiplying Eqs. (8) and (9) by t N , summing over N, and 
using the expressions in Eq. (10) for the corresponding coefficients of the terms in t and t 2 , 
we find that Z t , z[ A \ and z\ B ^ obey the following equations: 

z i = t + z\ A) + z\ B \ 



H) 

(l - z A t}z\ A) = z A t 2 Z t + z A t(l + t 
(l - z B t)z\ B) = z B t 2 Z t + z B t(l + t) . 



(12) 



Consequently, Z t is a rational function of t and is given explicitly by 



* = (13) 
C 2 (t) ' v ; 



where 



C 1 (t) = l + ZA + ZB -2t-t\ C 2 (t) = ^- l + ZA + ZB t + ? + t\ (14) 

Z A Z B Z A ZB Z A Z B 

Let ti, t 2 and *3 be the roots of C 2 {t) so that C 2 {t) = (t—t-\)(t—t 2 )(t—t 3 ). Expressing Eq. (13) 
in terms of elementary fractions and expanding the emerging factors l/(tj — t), j = 1,2,3, 
into Taylor series in powers of t/tj, Eq. (13) can be cast into the form 

where the coefficients are given explicitly by 

(i) = h + t 2 t 3 (2) = t 2 + ht 3 a{3) = t 3 + t x t 2 

(ti-t2)(ti-*3)' (t 2 - h)(t 2 - t 3 y (t 3 - h)(t 3 - 1 2 )' 

Equations (11) and (15) imply that the partition function of a chain with TV sites (Eq. (2)) 
is given explicitly by 

„ aW a 0) 

^ = ~jW + ~jW + TaT' ( 17 ) 

l l l 2 l 3 

and thus its behavior is completely determined by the properties of the roots of the cubic 
polynomial C 2 (t). Defining 

3 + 3z A + 3z B + z A z B 36 + 9z A + 9z B + 2z A z B r 
q= —. , r= , andX = — — , (18) 

54z A z B 9z A z B q 6 / 2 

one can show that Q = r 2 — q 3 < so that < X < 1 for all z AjB > 0, which implies that all 
three roots of C 2 (t) are real and given by [27] 

ii,3 = ±2^cos + i arcsin(X) j - ^ (19) 



i 2 = 2^/gsin Q arcsin(X)^ - i. 



(20) 



The roots of £ 2 {t) satisfy tit 2 *3 = — < 0, M2 + ht 3 + i 2 *3 = - - + — + ^ < 0, 

z A z B Z A Z B 

and since all three are real, it follows that one root is negative and the other two are positive. 
Using < arcsin(X) < tt/2, it is easy to prove that for any z A and z B one has t\ > t 2 > 
t 3 (thus t 3 is the negative solution) and 1*3 1 > t\. Consequently, due to Eq. (4), in the 
thermodynamic limit N — > oo the pressure is given by 

P = ~ In {l^ sin Q arcsin(X)) - ^ . (21) 
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We note that here three different roots appear since the model under study includes 
effectively three-site interactions, as shown by the indicator function <E>j in Eq. (3). The fact 
that as function of the chemical potentials neither pair of roots does intersect for any za and 
zb implies, as expected, that the one-dimensional model considered here does not exhibit a 
phase transition. More complicated models involving particles which require more than one 
empty lattice site for their adsorption, such as, for instance, the original ZGB model in which 
A species are dimers [1], can be treated within essentially the same approach as presented 
here. Such models would include effective four-site interactions (or more), and thus four or 
more different roots will emerge. 



4 Connection with classical spin S = 1 models 



Since, as we have noted, in the present model the hard-core interaction excludes the state 
(ni,rrii) = (1, 1) for a site, there is a natural connection between the present model and the 
three-state lattice gas model [28] or, equivalently, a special Blume-Emery-Griffiths (BEG) 
spin S = 1 model [29]. 

The mapping to the BEG model is accomplished as follows. We assign to each site of the 
one-dimensional lattice a three-state variable Oj, i = 1, . . . , N, such that 



+1, if the i-th site is occupied by an A particle, 
— 1, if the i-th site is occupied by a B particle, 
0, if the i-th site is empty. 

In terms of o^ the occupation numbers rij and rrii may be rewritten as 

ni = {<Ji + of )/2 , m, = (-<jj + of )/2. 



(22) 



(23) 



Periodic boundary conditions imply Ojv+i = <J\- Defining the nearest-neighbor (NN) coupling 
as (the limits corresponding to our particular model are indicated in parenthesis) 



J, 



-E 1 (- 0), 
-E 2 (- 0), 



for A-A neighbors, 
for B-B neighbors, 



+£?3 (— > oo), for A-B or B-A neighbors, 



(24) 



0, 



otherwise, 



one can write the Hamiltonian of this system as 

ft = [ n i n j(~ E i) + mimj(-E 2 ) + (riirrij + njm^E^] - ^(/x^n, + (iBm) (25) 



N 



<ij> 



i=l 
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where ^ means summation over all pairs of NN sites with each pair included only once. 



<ij> 



Replacing m, mi with the corresponding expressions (Eq. (23)) and collecting the terms, the 
Hamiltonian above may be rewritten as 



E\ + E2 + 2E3 s—^ E\ + E2 — 2E3 s—^ 22 E\ — E2 v-^ / 2 , 2\ 
n = 1 2^ 1 a i a o 1 + a 3°i ) 

(26) 



<ij> 



<ij> 



<ij> 



N N 
9 2^* 9 2^i' 



i=l i=l 

i.e., takes the form of a Hamiltonian for a spin 5 = 1 model [30] with the parameters 

E 1 + E 2 + 2E 3 E 1 + E 2 -2E 3 E x - E 2 „ ha - Hb , A MA + Mb 

= ; > ^ = ; , = ; , -n = z , and A = . 

4 4 ' 4 2 2 

Therefore, our model (E± = E2 = 0, E3 — > 00) is equivalent with a BEG model char- 
acterized by a rather unusual set of interaction parameters: C = 0, a bilinear exchange 
constant J = £3/2 — ► 00, and a biquadratic exchange constant K = —E3/2 — > —00 such that 
J/if = -1. 

The description in terms of the spin variables <7j opens an alternative way of computing 
the partition function in Eq. (2) by using the standard transfer-matrix method [31]. Although 
it is possible to set up the transfer matrix using directly the Hamiltonian given in Eq. (26), 
it turns out to be more convenient to define local fields //(<Tj) as 



-Ha, if <Ti = 1, 
+fJtB, if (Tj = -1, 
0, if <n = 0, 



(27) 



such that the partition function for the system is written as 



4 BEG) =£exp 
Wi} 



N 



i=l 



(28) 



or, equivalently, as the trace of the product of transfer matrices 



N 



(29) 



i=i 



where the transfer matrix V^j+i is given by 



= exp -/3Jj,i+icrjCTj + i - p(ji((7i)(Ti + ^(o- i+ i)cr m )/2 



(30) 
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In the thermodynamic limit the expressions for the pressure given by the partition functions 
in Eqs. (2) and (28) become identical if E\ = Ei = 0, and £3 — > 00. For these values of the 
parameters, the transfer matrix V£ j+i becomes 



Vn = V 



i,i+l 



ZA 
~A 





Z A 
1 





^z~B 

ZB 



\ 



(31) 



The eigenvalues Xj = 1/ tj of Vo are given by det(Vo - A/) = 0, which are the inverse of the 
zeroes of £2^) = 0. 



5 Mean density and compressibility 

5.1 Thermodynamics 

Since the model and the partition sum are symmetric with respect to interchanging A and 
B particles, it is sufficient to analyze the behavior of only one of the two densities, e.g., tia; 
the corresponding results for ng, and similarly for the compressibilities, follow by simply 
switching the indices A and B. 

Using Eq. (21), the mean density ua and the compressibility >ca of the phase formed by 
the adsorbed A particles are given as 

dP z A dn A , QoA 

n A = pz A -^ — , >c A = — ^ — • (32) 
Oza n A oza 

The behavior of ua as a function of za and zb, respectively, is shown in Fig. (2) for several 
values of zb and za, respectively. 

As can be seen in Fig. 2, the density ua is very small for za < z A (where (z A ,n A ) is 
the turning point of ua{za), occuring for zb > 1), raises sharply to large values for za > z* A 
and, independent of the fixed value for zb, approaches asymptotically ha(za — ► 00; zb) — ► 1- 
The functions z a (zb) and n A (zB) depicted in Fig. 3 show that the value z A is an almost 
linear function of zb, satisfying z* a {zb) ~ zb, while the density n* A = nA(z A (zs), zb) raises 
sharply with increasing zb and saturates at the asymptotic value ti a (zb — > 00) — > 0.5. The 
width Aza of the sharp increase in ua depends weakly on zb and at fixed zb the function 
tia(za'i zb) converges slowly toward the step function Q(za — zb) with increasing zb- In this 
limit the system exhibits two generic states: an ^4-depleted one in a "C 1, for za < z* A ) and an 
^-saturated one (ha ^ 1, for za > z A ), for which the density ua depends only very weakly 
on za, separated by a narrow transition region centered at za = z A . Similar conclusions can 
be drawn regarding the behavior of the mean density ua as a function of zb- Comparing 
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20 40 20 40 

Figure 2: Mean density ua as a function of za (a) and zb (b), respectively. From left to right, the 
curves correspond to (a) zb — 1, 5, 10, and 20 and to (b) za — 1, 5, 10, and 20, respectively. The 
dots indicate the turning points (z* A ,n* A ) and (z B ,n* A ). For comparison, the insets show results for 
zb = 10 (a) and for za = 10 (b), respectively, together with corresponding results for the case of no 
AB reaction (classical Langmuir, Eq. (5), dashed line) and for the approximation given by Eq. (6) 
(hard-core exclusion only, dotted line). 




these results with the ones corresponding to Langmuir adsorption (Eqs. (5) and (6)) it is 
evident that in general the reaction AB leads to a qualitatively different behavior of the mean 
density ua (see the insets in Fig. 2). We note that, as expected, at fixed za (inset Fig. 2(b)) 
the behavior for zb <C 1 is well described by the classical Langmuir non-interacting particles 
model but, rather counter-intuitively, the behavior at large za and fixed zb (inset Fig 2(a)) is 
also similar to that of a non-interacting case and not similar to the corresponding case with 
hard-core interaction. These findings will be addressed in the discussion of the asymptotic 
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behavior below. 

In Fig. 4 we show the compressibility of the phase formed by the adsorbed A particles 
as a function of the activities za and zb, respectively, for the same values as in Fig. 2. The 




20 40 20 40 



Z A 

Figure 4: Logarithm of compressibility ka as a function of za (a) and Zb (b), respectively. From 
left to right, the curves correspond to (a) zb = 1, 5, 10, and 20 and (b) za = 1, 5, 10, and 20, 
respectively. The squares denote the position of the turning points of ha{za) (Fig. 2(a)) and ha{zb) 
(Fig. 2(b)), respectively. 



behavior of the compressibility is consistent with that of the mean density, i.e., for example 
Fig. 2(a) shows that the low A-density phase at za < z* A has a large and slowly decreasing 
compressibility, while the high yl-density phase at za > z* A has a very small compressibility. 

As indicated by Eqs. (21) and (32), the functions ua(za,zb) and xa{za,zb) exhibit 
complicated dependences on za, zb, and thus an analytical study of the behavior of the 
system for all the values of the parameters is difficult. However, as we have mentioned, for 
fixed values of one of the activities, for example zb, there are three particular cases as a 
function of the other one, i.e., za' za <C z a , za 3> z A , and za — z* A . Below we derive and 
discuss the asymptotic behavior in these limits, starting from corresponding expansions of 
ua{za,z b ) or, if necessary, of t 2 (zA, z B ) (Eq. (20)). 

5.2 Asymptotic expansions 

We first consider the case when either one or both activities are small, i.e., the corresponding 
vapor pressure is low, the temperature is high, or the barrier against desorption is small. In 
the case when za "C 1, while zb is fixed, we find from Eqs. (21) and (32): 

1 (2z% + 6z B ~ 1) o / , N 

UA = 7V—^ ZA + TTT ,6 Z A + O (4) , (33) 



11 



and thus the compressibility is 

_ (l + z B f (4 + 64 + 11^ + 9) z B ^ , . . 

k a — z A + O [z A ) . (34) 

ZA (1 + Z B ) 

We note that in Eq. (33) the coefficient of the term linear in z A is smaller by a factor (1 + zb) 3 
than the corresponding one in the case of non-interacting Langmuir adsorption (Eq. (5)) and 
by a factor (1 + zb)~ 2 compared with hard-core Langmuir adsorption (Eq. (6)). Consequently, 
the compressibility is increased by similar factors. This reduction of n A is significant for large 
zb and shows that in the presence of many adsorbed B particles most of the adsorbed A 
particles do react with B and leave the chain. In the opposite case zb <C 1 the first-order 
term is the same as for the non-interacting Langmuir adsorption (Eq. (5)) since this limit 
corresponds to very low densities of both A and B particles and, consequently, very unlikely 
reaction events. 

Second, within the limit zb <C 1, while z A is fixed, we obtain 

ZA 3z A (44 + 13z A - 6)jA 2 , „ ( 3 \ /okn 

n A = — — ~ — 74 z B — -f^ z B + O (z B ) (35) 

l + z A (l+z A ) (l + z A ) 

and 

*A = j-+ 3( ^ +1 U + 0(4)- (36) 
ZA Z A (1 + ZA) 

As expected, the first term in the series for both ua and xa is the trivial Langmuir adsorption 
model result (Eq. (5); see also the inset in Fig. 2(b)), because the very low density of adsorbed 
B particles leads to a very small probability for a reaction AB. 
Third, for za = zb = z <C 1 we find 

n A = n B = z-Az 2 + 19z 3 + 0(z 4 ) (37) 

and 

x A = - + 3 + 3z - 8z 2 + O (z 3 ) , (38) 

which shows that at small activities the value of the density ua (or n B ) at the crossover 
za = zb increases linearly with the activity. 

Next we turn to the case when either one or both of activities are large, which can be 
realized in systems with low temperature, suppressed desorption, or at high pressures of the 
corresponding vapor phase. Most of the previous theoretical work has been focused on this 
limit [23]. Since in this case the analysis is somewhat more complicated compared to the 
situations with small activities, it is advantageous to start from the asymptotic behavior of 
*2 given in Eq. (20). 
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First, we consider the case za 3> 1 and z B fixed. Using the identity 

sin - arcsm JL - NO /0 = -\ , (39) 

V3 V 2(3 + z B )3/2 JJ 2y 3 + z B K J 

one finds for the Laurent series of £2 

z A z A z A z A \z A j 



This implies for the pressure 



1 1 , {l + 3z B ) 



^-* M+ «-^ + »ir +0 w (41) 

and, hence, for the mean density of the A particles 

»^-h+y^ +o (?Y (42) 

ZA Z A Z A \Z A J 

Therefore, in this limit the dependence on zb appears only in the third-order term and thus 
it is very weak. This explains the confluence of the density curves ha in Fig. 2(a) in the range 
of large za- Here we also note that regardless of the value of zb (provided that zb <C za) 
the mean density of the A particles tends to unity, i.e., the one-dimensional chain becomes 
saturated with A particles. Moreover, the first three terms in the expansion in Eq. (42) 
coincide with the corresponding terms in the expansion of the mean particle density for large 
za in the non- interacting Langmuir adsorption model (Eq. (5)). However, only the first 
term agrees with the similar expansion of the hard-core Langmuir adsorption (Eq. (6)). This 
confirms the numerical results shown in the inset in Fig. 2(a) and suggests that the reaction 
term exactly cancels the contribution of the hard-core interaction up to the third order in 
1/za- Intuitively, this can be understood by noting that for za S> zb the system tends to 
a state with a very low density n B , and thus the hard-core constraint that a site cannot be 
occupied simultaneously by an A and a B particle is effectively irrelevant, but the constraint 
that a site cannot be occupied by two A particles becomes very important. In this limit we 
obtain for the compressibility 

«4 + Tf+<>fi). <«) 



ZA Z A 

This confirms that in leading order the adsorption is given by the Langmuir-type result for 
the non-interacting adsorption model with corrections occurring only in third order. 
In the opposite limit > 1 and za fixed, similar calculations lead to 

t2 = ±-^-(l±fd> + ( 1+ ^- 4) +°(\), m 

zb z b z b z b z b \ z b / 
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which implies 




z A (3 - 2z A ) z A (6 - 6z A + 3z A ) z A 




) 



(45) 



and, respectively, 




+ O (z B ) . 



(46) 



Thus riA exhibits a very fast power-law decay implying a strong divergence of the compress- 
ibility xa for increasing zb- 

Lastly, for za = zb = z S> 1 we find 



5.3 Monte Carlo simulations 

The result in Eq. (47) tells that in this limit the sites are on average equally occupied by 
A and B particles. Since these findings give no insight into the actual spatial arrangement, 
i.e., the correlations, we have performed a simple Monte Carlo simulation. We start with 
an empty lattice of ./V sites and update the occupation numbers as follows. At any step, an 
empty (0) site is picked at random and an A (+1) oraB (-1) particle is deposited with equal 
probabilities (corresponding to equal activities za = Zb = z). If any of the two neighboring 
sites are occupied by a particle of opposite sign, the deposited particle and the neighboring one 
are removed, i.e., the occupation numbers of the sites are set to 0. If both of the neighboring 
sites are occupied by particles from the other species, one of them is randomly selected for 
desorption. Desorption of individual particles is disregarded since we are considering the limit 
z — > oo. Fig. 5 shows typical simulation results for the evolution of the system. The results 
clearly show the formation of domains consisting of mainly A or B particles such that only a 
small number, eventually just two in the limit N — > oo, of large domains emerges. Intuitively, 
the emergence of two domains is due to the fact the in the limit z — > oo the number of empty 
sites should be minimum since the ratio of the probability of a state with a greater number of 
empty sites and the probability of the state with minimum number of empty sites goes to zero 
as z — ► oo. Here we note that for a finite size system increasing the number of updates will 
eventually lead to only one domain of A or, equally probable, B, spanning the whole lattice. 
This is due to the fact that for z — > oo tiny deviations from za = zb, i.e., from exactly equal 
probabilities of depositing A or B, are sufficient to drive the system into one of the ua = or 
ha = 1 states (see also Fig. 2). These results show that in the limit za = zb — > oo, the system 
undergoes, in the thermodynamic limit N — ► oo, segregation in a state in which half of the 
sites belong to a domain of either A or B particles. Interestingly enough, similar behavior 




(47) 
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Figure 5: Occupation of a lattice of size N = 5000 after (a) 2.5 x 10 6 , (b) 5 x 10 6 , (c) 7.5 x 10 6 , and 
(d) 10 7 updates. Particles of type A and B correspond to +1 and -1, respectively. 



has been predicted for the steady-state behavior of diffusion- limited A + B — > reactions in 
systems with steady injection of the reactive species by external sources with equal intensities 
[32, 33] (for more details see also Ref.[21]). 

6 Conclusions. 

In this paper we have presented a model of monomer-monomer A + B — ► catalytic reactions 
on a one-dimensional chain in contact with a reservoir of A and B particles. The model 
assumes continuous exchange of A and B species between the chain and the vapor phase 
acting as a reservoir, and instantaneous reaction and desorption of neighboring A and B 
particles. We have calculated exactly the partition sum taking into account equilibrium 
fluctuations. From this we have obtained the pressure of the adsorbed particles, the mean 
density, and the compressibility. Altough in this one-dimensional model there is no phase- 
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transition, the system exhibits a rather non-trivial behavior. In particular, the mean density 
of the A particles changes rapidly from very small values when the activity zb is larger than 
the activity za to a state in which the lattice is occupied predominantly by A particles when 
za is larger than zb- In the case when the two activities are exactly equal and large the 
system undergoes segregation such that each of the two species clusters into a large domain 
and occupies half of sites of the chain. 
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